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The application of the nonperturbative renormalisation group approach to a sys- 
tem with two fermion species is studied. Assuming a simple ansatz for the effective 
action with effective bosons, describing pairing effects we derive a set of approximate 
flow equations for the effective coupling including boson and fermionic fluctuations. 
The case of two fermions with different masses but coinciding Fermi surfaces is con- 
sidered. The phase transition to a phase with broken symmetry is found at a critical 
value of the running scale. The large mass difference is found to disfavour the for- 
mation of pairs. The mean-field results are recovered if the effects of boson loops are 
omitted. While the boson fluctuation effects were found to be negligible for large 
values of ppa they become increasingly important with decreasing ppa thus making 
the mean field description less accurate. 
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The properties of asymmetric many fermion systems have recently attracted much at- 
tention (see, for example Ref. jl| and references therein) driven by the substantial advance 
in experimental studies of trapped fermionic atoms. This asymmetry can be provided by 
unequal masses, different densities and/or chemical potentials. Understanding the pairing 
mechanism in such settings would be of immense value for different many fermion systems 
from atomic physics to strongly interacting quark matter. The important theoretical issue 
to be resolved here is the nature of the ground state. Several competing states have been 
proposed so far. These include: LOFF Q phase, breached-pair (BP) superfluidity (or 
Sarma phase) and mixed phase |4] . Establishing the true ground state is still an open ques- 
tion. It was shown, for example, that LOFF and mixed phases are more stable then the 
Sarma phase in the systems of fermions with the mismatched Fermi surfaces and with both 



equal and different masses jjj, |4|, |5[ . All these studies, however, have been performed within 
the mean field approximation (MFA). In spite of the fact that in many cases MFA is quite 
reliable it is important to understand better the limits of applicability of MFA and work out 
the physical regimes where the MFA is too crude or even inadequate. The convenient way to 
estimate the corrections to MFA is provided by the nonperturbative renormalisation group 
(NRG) approach [f| which was successfully applied to the standard pairing problem with 



one type of fermions ^1 [icj| . The main element of NRG is the effective average action 
which is a generalisation of the standard effective action T, the generating functional of the 
1PI Green functions. The only difference between them is that Tk includes only quantum 
fluctuations with momenta larger then the infrared scale k. The evolution of the system as 
the function of the scale k is described by the nonperturbative flow equations. When k — > 
all fluctuations are included and full effective action is recovered. Similarly, at starting scale 
k = K no fluctuations are included so T k= x can be associated with the classical action 
S therefore provides an interpolation between the classical and full quantum effective 
actions. 

The dependence of from the infrared scale k is given by the nonperturbative renor- 
malisation group equation (NRGE) 

d k T = - t -Ti[(d k R)(T^-Rr 1 }. (1) 

Here r*- 2 -* is the second functional derivative of the effective action taken with respect to 
all types of field included in the action and R(q, k) is a regulator which should suppress 
the contributions of states with momenta less than or of the order of running scale k. To 
recover the full effective action we require R(q, k) to vanish as k — > whereas for q « k the 
regulator behaves as R(q, k) ~ k 2 . The above written equation is, in general, the functional 
equation. For a practical applications it needs to be converted to the system of partial or 
ordinary differential equations so that approximations and truncations are required. 

We consider a nonrelativistic many-body system at zero temperature with two types of 
the fermion species a and b interacting through a short-range attractive interaction and 
introduce a boson field describing the pair of interacting fermions. The ansatz for T takes 
the form 
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T[ip,^,(j),(j)\p,k} = / d A x 



- Z 9 [ \ ^b^atf ~ ^ ^l^Ft 



<t>{x)-u{<t>A ] ) 



(2) 



Here Mj is the mass of the fermion in vacuum and the factor l/2m with m = M a + M b in 
the boson kinetic term is chosen simply to make Z m dimensionless. The coupling Z g , the 
wave-function renormalisations factors and the kinetic-mass renormalisations factors 
Zm t M all run with on k, the scale of the regulator. Having in mind the future applications 
to the crossover from BCS to BEC (where chemical potential becomes negative) we also 
let the chemical potentials p a and /if, run, thus keeping the corresponding densities (and 
Fermi momenta pp,i) constant. The bosons are , in principle, coupled to the chemical 
potentials via a quadratic term in 0, but this can be absorbed into the potential by defining 
U = U — (fj,i + n-i)2i 0</>V- We expand this potential about its minimum, 0^0 = p , so that 
the coefficients Ui are defined at p = p , 



U (p) = M + ui{p - Po) + ^u 2 (p- po) 2 + - u 3 (p - p ) 3 + 



(3) 



where we have introduced p = (jy<j>. Similar expansion can be written for the renormalisation 
factors. The coefficients of the expansion run with the scale. The phase of the system is 
determined by the coefficient u±. We start evolution at high scale where the system is 
in the symmetric phase so that u\ > 0. When the running scale becomes comparable 
with the pairing scale (close to average Fermi-momentum) the system undergoes the phase 
transition to the phase with broken symmetry, energy gap etc. The point of the transition 
corresponds to the scale where U\ = 0. The bosonic excitations in the gapped phase are 
gap-less Goldstone bosons. Note, that in this phase the minimum of the potential will also 
run with the scale k so that the value p (k — > 0) determines the physical gap. 
The evolution equation takes the following general form 



ar = -- Tr 



(d k R B ) (rgi - Rb)- 1 ] + l -Tr [(d k R F ) (F% - Rp)-^ 



(4) 
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Here T^b^ff) * s ^ ne m &trix of the second functional derivatives of the effective action taken 
with respect to boson(fermion) fields included in the action and R B (R F ) is the boson 
(fermion) regulator which should suppress the contributions of states with momenta less 
than or of the order of running scale k. The boson regulator has the structure 

R B = R B diag(l,l). (5) 

The fermion regulator for both types of fermions has the structure 

Rf,i = sgn(ei(q) - Hi)R F:i (q, fi t , k)diag(l, -1) (6) 

Note that this regulator is positive for particle states above the Fermi surface and negative 
for the hole states below the Fermi surface. 

The evolution equations include running of chemical potentials, effective potential and all 
couplings (Zfj,, Z m , Z M>i , Z^j, Z g ). However, in this paper we allow to run only Z^, parameters 
in the effective potential (u's and po) and chemical potentials since this is the minimal set 
needed to include the effective boson dynamics. 

Calculating the second functional derivatives, taking the matrix trace and carrying out 
the pole integration in the loop integrals we get the evolution equation for U at constant 
chemical potentials 

d k U = -^-d k T = ~ f —JL F ' S [sgn(g - p^ a ) d k R F , a + sgn(g - p^ b ) d k R F , b ] 

V 4 2 J W Je 2 FiS + A 2 

1 f d 3 q E BR 

+ ^;Jw?7W^l dkRB - (7) 



Here 



and 



E s = (Ep :a + E FJb )/2, E A = (E F , a - E b , a )/2, (8) 



E B {q,k) = 1 ^q 2 + u l +U2{2^<j)-p )+RB{q,k), V B =u 2 ^<f), (9) 
E F ,i{<liP»,ii k ) = ^f.V 2 ~ + RF(q,k) sgn(g-p^), A 2 = g 2 ^(f). (10) 



and we have introduced p^i = y/2M i /j, i , the Fermi momentum corresponding to the (running) 
value of Hi. It is worth mentioning that poles in the fermion propagator occur at 



\* = -E A ±y/Efak? + &- (11) 



5 



At k = (Rf = 0) in the condensed phase, these become exactly the dispersion relations 
obtained in [3| where the possibility of having the gapless excitations has been discussed. 
The ordinary BCS spectrum can easily be recovered when the asymmetry of the system is 
vanishing (Ea — > 0). The first term in the evolution equation for the effective potential 
describes the evolution of the system related to the fermionic degrees of freedom whereas 
the second one takes into account the bosonic contribution. The mean field results can be 
recovered if the second term is omitted. In this case the equation for the effective potential 
can be integrated analytically. 

d 3 q 



U(p,p,k) = U(p,p,K) 



(27T) 



y/E s {q, kf + A 2 - y/E s (q, Kf + A 2 



At starting scale K the potential has the form 



U(p,p,K)=u (K)+ Ul (K)p. 
The renormalised value of U\(K) can be related to the scattering length. 



u-l(pf,K) 
9 2 



M 
2na 



(27T) 
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|£ s (g, 0,0,0) Es(.q,tia,Hb,K)\ 



;i2) 



(13) 



(14) 



Here M is the reduced mass and the dependence of Eg on the chemical potentials has been 
made explicit. 

Differentiating the effective potential with respect to p, setting the derivative to zero and 
taking the limit K — > oo, we arrive at the equation 

d 3 q 



M 1 
2vra + 2 



(27T) 



0. 



(15) 



Es{q, 0,0,0) ^E s (q,p a ,p b ,ky + A\ 

Taking the physical limit (k = 0) we obtain the gap equation identical to that derived in 
the mean field approximation Q, Q]. 

We now turn to the full set of the evolution equations which includes the effects of the 
bosonic fluctuations. In this paper we consider the case of two fermion species with the 
different masses and the same Fermi momenta. It implies that the chemical potentials are 
different. In this situation the Sarma phase does not exist and the system experiences the 
BCS pairing depending however on the mass asymmetry. The general case of the mismatched 
Fermi surfaces will be discussed in the subsequent publication. 

The derivation of the evolution equations was discussed in details in Ref. Q so that 
here we just mention the main points. Within the above described approximation (fixed 
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couplings Z m , Z M>i , Z^ ti , Z g ) all of these can be obtained from the evolution of the effective 
potential, for example 

1 d 2 



^ ' p=po 

where \i = p, a + p,b- Substituting the expansion for the effective potential on the left-hand 
side of the evolution equation leads to a set of ordinary differential equations for the running 
minimum p and coefficients u n . These equations have a generic form 



d k u n - u n+1 d k p 



dp' 



(d k u) 



(17) 



P=po 



One can see from this equation that some sort of closure approximation is needed as the 
equation for u n always include u n+ \ coefficient etc. In this paper we calculated u n> 2 in the 
MFA with the effective potential given by the Eq.(12). As already mentioned we follow the 
evolution of the chemical potential keeping density fixed. Defining the total derivative 



— - | dp0 9 
dk k dk dpo 



(18) 



and applying it to the ^ (or to |n) we obtain the following set of equations 

dp, d 



dpo 
dk 



d,U 



(19) 



P=po 



dk dp, 

where is the coefficient in the leading term of the expansion for Z^ similar to Eq.(3), and 

d k U\ 



du da 

h n — 

dk dk 



I p=po 





-u 2 


dpo 
dk 


du 2 


- u 3 


dp 


~dk 


dk 


dz^o 
dk 


- Z^i 


dpo 
dk 
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dp, 
dk 

dp, 
dk 



dp 2 



d k U 



P=PO 



P=PO 



1 , dp 

2 X dk 



1 d 2 



2 dfxdp 



(20) 
(21) 

(22) 

(23) 



P=P0 



where we have defined 



X 



d 3 U 



dp?dp 



Z(j,l 



1 d 3 U 



P=P0 



2 dp,d 2 p 



(24) 



P=P0 



These functions have also been calculated in the MFA. The set of evolution equations in 
symmetric phase can easily be recovered using the fact that chemical potential does not run 
in symmetric phase and that p = 0. 
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Let us now turn to the results. For simplicity we consider the case of the hypothetical 
"nuclear" matter with short range attractive interaction between two types of fermions, 
light and heavy, and study the behaviour of the energy gap as the function of the mass 
asymmetry We choose the Fermi momentum to be pp = 1.37/m -1 . One notes that the 
formalism is applicable to any type of a many-body system with two fermion species from 
quark matter to fermionic atoms so that the hypothetical asymmetrical "nuclear" matter is 
simply chosen as a study case. We assume that M a < M&, where M a is always the mass of 
the physical nucleon. 

In this paper we use a sharp cutoff function chosen in the form which makes the loop 
integration as simple as possible 

k 2 

RF,i = ^ M ~M k + P^) 2 - <? 2 )%am + k-q) + ((k + p^) 2 + q 2 - 2pl-)6(q - p^ + k)] , 

(25) 

and similarly for the boson regulator 

R B = ^(k 2 -q 2 )9(k-q). (26) 
iere 6{x) is the standard step- function. This type of boson regulator was also used in Ref. 



ll| (see also Ref. 

The use of a sharp cutoffs can be potentially dangerous as it may generate the artificial 
singularities when calculating the flow of the renormalisation constants (Z's) but seem to 
be harmless when all the evolution parameters are related to the effective potential RG flow 
as is the case here. 

As we can see the fermion sharp cutoff consists of two terms which result in modification of 
the particle and hole propagators respectively. The hole term is further modified to suppress 
the contribution from the surface terms, which may bring in the dangerous dependence of 
the regulator on the cutoff scale even at the vanishingly small k. We found that the value of 
the gap practically does not depend on the starting point provided M a ^ « K. As expected, 
the system undergoes the phase transition to the gapped phase at some critical scale which 
depends on the value assumed for the parameter ppa where a is the scattering length in 
vacuum. One notes that the critical scale does not depend on the mass asymmetry. 

First we consider the case of the unitary limit where the scattering length a = — oo. The 
results of our calculations for the gap are shown on Fig. 1. 
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FIG. 1: Evolution of the gap in the MF approach (dashed curve) and with boson loops (solid 
curve) in the unitary regime a = — oo as a function of a mass asymmetry. 

We see from this figure that increasing mass asymmetry leads to a decreasing gap that 
seems to be a natural result. However, the effect of the boson loops is found to be small. 
We found essentially no effect in symmetric phase, 2 — 4% corrections for the value of the 
gap in the broken phase and even smaller corrections for the chemical potential so that one 
can conclude that the MF approach indeed provides the reliable description in the unitary 
limit for both small and large mass asymmetries. It is worth mentioning that the boson 
contributions are more important for the evolution of u<i where they drive «2 to zero as 
k — > making the effective potential convex in agreement with the general expectations. 
This tendency retains in the unitary regime regardless of the mass asymmetry. 

We have also considered the behaviour of the gap as the function of the parameter p F a 
for the cases of the zero asymmetry M a = M b and the maximal asymmetry M b = lOM a . 
The results are shown on Fig. 2. 

One can see from Fig. 2 that in the case of zero (or small) asymmetry the corrections 
stemming from boson loops are small at all values of the parameter p F a considered here 
(down to p F a = 0.94). On the contrary, when M b = 10M a these corrections, being rather 
small at p F a > 2 becomes significant (~ 30%) when the value of p F a decreases down to 
p F a ~ 1. We found that at p F a ~ 1 the effect of boson fluctuations becomes non negligible, 
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FIG. 2: Evolution of the gap as a function of the parameter ppa. The upper pair of the curves 
corresponds to the calculations with no asymmetry in the MF approach (dashed curve) and with 
boson loops (solid curve) and the lower pair of the curves describes the results of calculations with 
the maximal asymmetry when M& = lOM a 

~ 10% already for M b = 5M a . One can therefore conclude that the regime of large mass 
asymmetries, which starts approximately at Mb > 5M a , moderate scattering length and/or 
the Fermi momenta is the one where the MF description becomes less accurate so that the 
calculations going beyond the MFA are needed. One might expect that the deviation from 
the mean field results could even be stronger in a general case of a large mass asymmetry 
and the mismatched Fermi surfaces but the detailed conclusion can only be drawn after the 
actual calculations are performed. 

We were not able to follow the evolution of the system at small gap (or small p F a) 
because of the non-analyticity of the effective action in this case which means that the 
power expansion of the effective potential around the minimum is no longer reliable. To find 
the evolution at small gap the partial differential equation for the effective potential should 
probably be solved. 

To summarise, we have studied the pairing effect for the asymmetric fermion matter with 
two fermion species as a function of fermion mass asymmetry. We found that regardless of 
the size of the fermion mass asymmetry the boson loop corrections are small at large enough 
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values of ppa so that the MFA provides a consistent description of the pairing effect in 
this case. However, when ppa ~ 1 these corrections become significant at large asymmetries 
(Mb > 5M a ) making the MFA inadequate. In this case it seems to be necessary to go beyond 
the mean field description. 

There are several ways where this approach can further be developed. The next natural 
step would be to consider the case of the mismatched Fermi surfaces taking into account the 
possibility of formation of Sarma, mixed and/or LOFF phases and exploring the importance 
of the boson loop for the stability of those phases and applying the approach to the real 
physical systems like fermionic atoms, for example. Work in this direction is in progress. 
The other important extension of this approach would be to include running of all couplings 
of the effective action and use different type of cut-off function, preferably the smooth one. 
The three body force effects Q], when the correlated pair interact with the unpaired fermion 
may also turn out important, especially for non-dilute systems. 

The author thanks Mike Birse, Niels Walet and Judith McGovern for numerous very 
helpful discussions. 
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